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ABSTRACT 


The transient behavior of combat logistics support systems is analyzed. 
Combat availability is defined as the number of active combatant platforms being 
supported by a single fault diagnosis and repair facility. Heavy traffic conditions 
inherent to intense combat periods allow the use of diffusion approximation models, 
which provide speedy solutions used to compare adaptive scheduling policies to a 
standard First-Come, First-Serve policy. The adequacy of these models is investigated 
and numerical solutions are compared to simulation results. The case in which failed 
modules require a degree of support that is beyond the capability of local 
maintenance is also investigated for both pre- and post-local-repair relocation to 
distant repair. The use of cannibalization in short-term situations is shown to have 
a dramatic effect in terms of combat availability. A preliminary model for a non- 
cannibalization policy is also developed. Optimization models for choosing spare 
parts allocation within budget constraints, or for achieving required operational 


availability with minimum cost are described. 
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I. INTRODUCTION 


A. THE NATURE OF THE PROBLEM 

At important times in military operations, especially during intense conventional 
combat periods, the problem of maintaining the readiness and operational availability 
of combatant units is of great relevance. Usually the available repair facility is 
restricted in terms of diagnosis and/or maintenance capability, since it also 
constitutes part of the deployed group. Thus a detachment unit is expected to operate 
as successfully as possible for a certain period of high activity, during which relatively 
many equipment failure events are likely to occur. The purpose of this thesis is to 
show how the desired availability can be enhanced. 

As an example, one can imagine a Carrier-based Air Unit the aircraft of which 
experience diverse failures of different mission-essential modules; the failures can be 
either total or partial. A total failure indicates that the original parts inventory may 
be depleted permanently with the removal of these components from service, while 
partial depletion signifies that these failed modules are repairable: they can either 
be immediately served or else join a queue for future service. Also, there is the case 
in which a failed component arriving at the initial diagnosis station, even if not 
rendered completely useless, is indicated to require a degree of support that is 
beyond the capability of local maintenance ("BCM"). The alternative may be to send 


this failed module to other maintenance levels, with consequent transit and repair 


times uncertainty being introduced. Because of the delays involved it is clearly 
important to identify as BCM only those failed components that truly need the 
distant Service. 

A basic and important operational problem can be summarized as follows: how 
to determine a good schedule for repairing failed modules in order to maximize 
combat availability, subject to various resource restrictions? Usually there is no 
fundamental reason for servicing failed units in the order in which they fail ("First- 
come,First-serve” policy), although it is natural and superficially "democratic" or "fair" 
to do so. Other service policies, allowing for queue length influence, may very well 
grant increased system availability as a function of time, defined as the total number 
of operational units "up" at any time f. 

This thesis develops and exercises various mathematical models for evaluating 
scheduling and spares stockage rules in a transient dynamic combat environment. It 
is very important that adaptive combat logistics support models have the ability to 
produce a description of transient behavior, since a steady-state situation may never 
be achieved. It is also desirable to get a solution for both the mean number of 
components in the system at each time ¢ and the variances as well, since this 
knowledge will permit the use of various measures of effectiveness which are closer 
in meaning to standard definitions of availability, i.e., the probability that the number 
of available units, e.g. carrier-based aircraft, exceed a pre-specified value at each 


time f. 


In this study operational units will be referred to as aircraft, and modules may 
be considered as being major avionics components, distinct in each aircraft, but 
assumed to be essential to the operation. It is important not to let the generality of 
the problem be obscured by these considerations. 

Stochastic models for these distributed logistic systems have been derived and 
verified via simulation for situations such as one server at the repair facility, and 
Markovian failure and repair times. These analytical models may be employed as 
tools to analyze the effect of different service disciplines, instead of the exclusive use 
of time-consuming and large computer-intensive simulation techniques, such as 
DYNAMETRIC, a package used at RAND Corporation. It is emphasized that the 
development of the present analytical modeling methodology is in its infancy, and 


that DYNAMETRIC remains a standard valuable tool. 


B. DEFINITIONS 

Mathematical definitions and formulations follow. 

Let the index 7 identify each of the different types of components to be 
considered; 7 = 1,...,/. The Aircraft Unit, together with its base of operation, deploys 
with K; components of type 7, considering components effectively installed in the 
aircraft plus any spares to be kept on base, e.g. in the Aircraft-Carrier local parts 
inventories, to replace those completely lost through attrition or to keep the aircraft 


operational while failed components undergo repair or wait for service. K; may 


actually change over the period of combat as some items permanently fail but others 
are added. 

A local maintenance shop is available to provide service/repair for failed 
components. It is assumed that components have independent exponential failures 
at constant rate 4; and that times to repair are also independent exponential with 
rate v,; our methods will actually accommodate time-dependent A; and v;. The single 
repairman at the repair shop "sees" arrival rates of failed components which are 
equal to individual failure rates multiplied by the total number of items operating at 
each particular time. This number is, obviously, the number of aircraft actually 
operational at that time. At any time, ¢ let V(t) be the number of components of type 
i waiting in queue or being serviced. 

Let A, be the number of aircraft initially deployed, and let A,(t) denote the 


number of operational (i.e., available) aircraft at time ft. A(t) is defined by 
Aft) = min {4,,K,-N,(0,K,-N,(0,..,K,-N«O) . (1.1) 


As defined earlier, the time-dependent failure rate of each type of component, 
as seen by the repairman, is equal to \,A,{t), for i=1,2,...,7. At each time, f there are 
K, - N,(t) components of type 7 available for use. The definition of A,(t) in (1.1) 
follows from the fact that the least available item determines the total availability of 
operational aircraft, if less than A,. Note that in this formulation A, can be a function 


of time if attrition due to combat is considered. 


Equation (1.1) assumes that only one item/module of each type is installed in 
each of the A, aircraft, and that cannibalization is allowed, i.e., an inactive (due to 
failure in component type /) aircraft may be removed from its working components 
in order to permit another unit, which suffered a failure of a different component, 
to return to active status. 

Another assumption that will be needed later is that the system is in heavy 
traffic, 1.€., 


I 
min {4,K, ,K,,.K)}x }o — > 1. (1.2) 


A; 
iz1 V; 

This assumption appears to be very realistic for most deployments, specially in 
a real combat situation. 

Considering now the repairman side of the problem, it seems clear that it is 
costly, in terms of time, to switch from the component being currently repaired to a 
new one before service termination and, consequently, the service strategies to be 
considered will determine the next component to be serviced (to completion) at the 
moment of the previous service completion. Of course if a server incumbent is not 
finished for an extraordinarily long time it may well be desirable to interrupt its 
Service and substitute another now more vital item. Rules for such a substitution are 


not evaluated in this thesis. It is clear that under some circumstances such procedures 


can be useful adaptations. 


C. GOALS AND SCOPE 

The present work attempts to exploit stochastic models developed for similar 
systems, as well as to perform initial sensitivity analysis with respect to certain tuning 
parameters. 

Monte Carlo simulations are employed as validation tools, and particularly 
extreme numerical examples are investigated so as to build understanding of the 
capabilities of this approach. 

Also, aS mentioned earlier, the original models are augmented in order to be 
able to consider the possibility of module failures which are "beyond the capability 
of maintenance". The objective is to show that, with small modifications, these 
models are capable of providing adequate approximations for the somewhat more 
realistic "BCM" problem. 

The effect of not using the cannibalization policy is investigated via a simulation 
model, and a modification to the analytic model is proposed to account for this 
change in strategy. 

Finally, an optimization model is proposed, allowing for considerations of 
limited resources and for the need of quick allocation of such resources to spares for 
the various modules, accounting for the effect of following an adaptive repair 


scheduling policy, so as to maximize overall system availability. 


Il. REVIEW OF EXISTING MODELS - BACKGROUND 


Models have been developed to analyze the transient behavior of similar 
systems with respect to the selection of a particular service policy. Initially, a natural 
choice for a Measure of Effectiveness ("MOE") is given by the expected number of 
aircraft available at time ft, E[A,{f)], (1.1) or, nearly equivalently, by the set of 


expected values of available modules of each type, E[K; - N;(t)], 7=1,2....,/. 


A. SIMULATION 

Several maintenance policies have been compared in a previous work by Latta 
(Ref. 1]. Simulation was employed to study six different repair policies going from the 
original First-Come, First-Serve ("FC,FS") scheme to what can be called Least 
Available Item Next ("LAIN") policy. This last scheduling strategy determines service 
priority based on the current availability of all types of components, and in which the 
repairman, after scanning inventory levels (if any), re-orders components in the 
service queue to favor those with the lowest operational stock level. If all initial 
stocks are completely depleted at any time ¢, the component type with lowest value 
of K; - N,(t) is chosen to be serviced next. 

Latta demonstrated that the LAIN policy yields considerable improvement in 


terms of mean number of available aircraft over all other considered policies. 


A clear disadvantage of the simulation technique for exploring our model 
implications is the intensive requirement for computer usage inherent in the 
straightforward Monte-Carlo simulation approach. It seems prohibitively costly, if not 
entirely impractical, to employ these models in a optimization context, for example, 
or in cases where imprecisely determined parameters for failure and repair times 
indicate the need for some kind of dynamic assessment of uncertainty in these 
parameters, such as by "bootstrapping". Real problems involve aspects of both of the 


issues; this paper has something to say about each. 


B. STOCHASTIC MODELS 


1. Definitions 
Diffusion processes are particular cases of Gaussian processes with 
continuous sample functions, and were originally used to model physical phenomena, 
e.g. the motion of small particles of gas. The Brownian Motion (or Wiener process) 
{W, , t>0} is the most elementary example of a diffusion process, having zero drift 
(or infinitesimal mean), and diffusion parameter (or infinitesimal variance) 
independent of W. It is the continuous time version of a random walk. 


A diffusion process {X, , t>0} satisfying the stochastic differential equation 
dX, ==pX, + odW, Gens, = Cateye (2.1) 


is called an Ornstein-Uhlenbeck process. This stochastic process was first used to 


directly model the velocity of a particle subject to elastic forces. In this case, the drift 


parameter reflected a restoring force proportional to the distance, and directed 
towards the origin. 

For a complete treatment of diffusion processes see, for example, Karlin 
and Taylor [Ref. 2]. Arnold [Ref. 3] examines the multivariate stochastic differential 
equations for the Ornstein-Uhlenbeck process, and gives general results for the mean 


and covariance of X, in (2.1). 


2. Air Unit Detachment Model and Related Work 

Several papers have discussed the use of diffusion processes as an 
approximation for job-server systems under heavy traffic. See, for example, Gaver 
and Jacobs [Ref. 4], Gaver and Lehoczky [Ref. 5 and 6] and Iglehart [Ref. 7]. 

Gaver and Jacobs [Ref. 4] developed diffusion approximation models for 
computer systems with processor-shared service disciplines. Gaver and Lehoczky 
[Ref. 6] studied a repairman problem with two types of repair. Pilnick [Ref. 8] 
extended these models to account for multiple types of queues and service priority 
proportional to a function of queue length and modeled the Air Unit detachment 
problem directly. Gaver, Isaacson and Pilnick [Ref. 9] exploit these models and 


presented various applications. The results are summarized here. 


a. Processor-sharing Adaptation 
If we define g,(N(t)), where N(t) is the vector with components N,,(t), 


as the proportion of time jobs of type i are served by the processor, we can view 


q,(N(t)) as the probability that the processor (server) will select job 7 for service just 
after each time slice, when a departure takes place in that time slice of length dt. 

In order to derive mathematical models for the actual Air Unit 
problem, g,(N(t)) is used to represent the probability that, after service completion 
at time t, a module of type 7 is selected for being serviced next. In this case, it turns 
out that g,(N(t)) is, in fact, a function of N(t) and A,{Z), 1.e., it depends on the queue 
sizes and the operational (combat) availability. In fact, g,(N(t)) is a decision variable, 
subject to determination by the scheduler. For notational simplicity, denote g,(N(t)) 
by q;(t). 

6b. Diffusion Approximation - The Air Unit Detachment Problem 
Pilnick [Ref. 8:p.143] showed how the following system of stochastic 


differential equations can be obtained: 


dN(t) = A;Ay(t) dt - v,p,(t) dt 





| a0 
+ |A Ault) + vp) (1 + 2p,(O{v, Y EO eee ie) ae 
J 


Vj 

where {W.(t), t>0} are independent standard Wiener processes, N,(t), 1=1,2,...,J are 
continuous approximations of the actual discrete processes, and 

Ww. 

pt) = C—qt) , 

B; 
where C is a normalization constant depending on the specific form of p(t), and 
w,;>0 denote arbitrary scalars representing weights for items of type i; see Chapter 


III-A below . 
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Assume that the system is large, for example let a=A,+ 2K, , A,=a.a,, 
K,;=a.a;, and #;=a.v; are constants. Define B,(t)=N,(t)/a, where B(t) approaches a 
deterministic function of t representing the scaled mean of N,(t), as a. 


Consider the normalization/approximation 


N(t) - a B,{) 


X(t) « (2.3) 


a 


(1) Expected Values. If we define dN,(t) as the increment N,(t+dt)- 
Nt) and express it in terms of transformations similar to those in (2.3), it can be 
shown (see Pilnick [Ref. 8]) that the following system of / ordinary differential 


equations can be derived for 8,(t), 7=1,2....,/: 


dB ,(t) 


Ti = i, ay(t) - ppt) (2.4) 





where a,{t) is a scaled version of A,{t) in (1.1), Le. 
a,(t) = min{a,,a, - B,(),...,«,-B,@} (2.5) 


and q,(t) are smooth functions chosen to be a representation of service policy. 


Different definitions for q,(t) will be described in Chapter III. 


(2) Vanances. A scaled variance-covariance matrix of N,(t), 2(t) can 
be also defined through a system of /xJ ODE’s. Following the notation in Pilnick 


(Ref. 8, pp. 85,147], let 
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Oeil) Oral) ene 4D) 
T(t) « - — a ier 


Loy (t) ap(t) ... on(t)| 


be the variance-covariance matrix of N(t) scaled by (2.3). The equations for the 


elements of X(t) are 





I 
doa;;(t) _ Bone) + 25 H,6,(t) (2.6) 
dt = 
and, fori #j], 
do;. I 
se - S [ H;,(¢) G;,(t) + A,,(t) o,,(t) ] : (2.7) 
k=] 


Here, H(t) is the /x/ matrix 


H(t) ifj + s(t) 


H.(t) = 
yj HG) =e ie 


s(t) = argmin { 0, - Bt), ..., % - B,(t) } , with 
|: AC) an at CY C9 2.8 
p Es pBay” )(1 -p,)) (2.8) 


and, for # J, 
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Hoek eae eee y(t) p(t) (2.9 
NST YS aie . 


The scalars €; , p and y are determined by the choice of the service policy, as we will 


see in Chapter III. 


B?(t) in (2.6) is the JxJ diagonal matrix with elements 


pA _ 


B?,(t) = A,ay(t) + ppt) (1 + 2p, (pn PE - 13) . 
J 


J 
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II. ADAPTIVE SUPPORT LOGISTICS 


Chapter II introduced the general analytic expressions for the mean and the 
variance-covariance function of an approximation to N(¢), valid for a large system 
("large" suitably defined) in heavy traffic. In this chapter, adaptive scheduling policies 
are studied, and compared to a standard FCFS policy. The accuracy of the diffusion 
approximation is verified in various numerical examples by comparison with results 
of Monte Carlo simulations. The behavior of the analytic models with respect to 
certain tuning parameters is demonstrated in cases showing wide variation of repair 
and failure rates. Under suitable conditions the analytical approximation is in 
excellent agreement with simulation, and is conducted in a fraction of the computer 
time required by present simulations; see below. 

Most of the numerical solutions of the differential equations were obtained on 
an IBM 3033/4381 computer at the Naval Postgraduate School using the IMSL 
(Release 10) package’s ODE solvers with either Adam’s method or Gear’s stiff 
method (backward differentiation formula). All codes were written in such a way that 
immediate translation to PC-FORTRAN is possible. 

All simulations were carried out on the IBM 3033/4381 mainframe, using the 
LLRANDOMII random number generator package [Ref. 10]. Certain relevant details 


of the simulations are provided when necessary in the next sections. 
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It should be mentioned that, for 1000 replications, the time spent in the 
simulation models was, on average, on the order of eight minutes, while the 
numerical solutions of the differential equations usually took approximately one 


second for the examples under study. 


A. MODELING ADAPTIVE POLICIES 

As was mentioned in Chapter II, in the diffusion approximation model is 
essential to characterize the availability and queueing behavior induced by a 
scheduling policy represented by the choice of q,(t). Pilnick [Ref. 8] describes fairly 
general expressions for q,(t) that can be useful for modeling several service policies. 
Put simply, let 


pt) = Ww; q(t) /p,; 


ee os (3.1) 
af 


where q,(t) is the "probability" (a relative measure) that an item of type 7 is selected 
for service next, if the previous service is completed at t. The following fairly general 
form of q,(t) incorporates a plausible definition for the parameters p, y and &; in 


(2.8) and (2.9): 
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(&,+ p BM) ) 


1.0)= eee 
(+p ROY 
J 


(3.2) 


The actual numerical values assigned to the parameters will depend upon the 
selection of the scheduling policy. For our present purpose, two forms for q,(t) are 


investigated. 
1. First-Come, First-Served (FC,FS) 
A technique to model this common policy is to suppose that the probability 
of choosing a module of type 7 for being serviced next is proportional to N,(t). In this 
case we must put p=1, y=1 and €;=0 and (3.2) becomes 


B (2) 


At) = ’ 
WO EG (3,3) 
J 





2. Least Available Item Next (LAIN) 
Gaver, Isaacson and Pilnick [Ref. 9] devised the term “anti-availability" 
of a module i to characterize the quantity (K,-N,(t))". If q,(t) is defined as 


[a,- BC]? 


Oe pea 
> [a;-B,(t)]? 
J 


(3.4) 


where p is a large integer (e.g., 10 - 40), then (3.4) is an attempt to emulate a 
deterministic choice of the module with largest anti-availability. This present form 


of q,(t) implies that we must have c=-1, y=-p and €,;=a, in (3.2). 
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a. Sensitivity Analysis 

It turns out that the selection of the tuning parameter p in (3.4) 
influences the accuracy of the diffusion approximation, depending upon the range of 
values for failure and repair rates. 

In this section a hypothetical cases are examined: the first (Case 1) 
exhibits comparable values for the rates and it is demonstrated that low values of p 
yield good precision; the second case (Case 2) illustrates the fact that a larger value 
of p must be used when the rates vary considerably. For both cases in this section, 
A, = 50 and / = 10, and no spares are provided. It is assumed that we are interested 
in evaluating combat availability for missions of duration up to T=100 (e.g., days). 

Simulation results are used to establish the basis for comparison. 
Here, as well as in all numerical examples for the LAIN policy in this thesis, the 
simulation model chooses the actual least-available module for service next as it 
moves along the sample path in each replication. The analytical method represents 
a probabilistic selection, but one that with near certainty picks the least available 


item for repair. 


IF 


(1) Case 1. Table 3.1 displays the input data for the first case. 


Failure and repair rates are given with respect to a "standard" time unit (e.g., days). 


Table 3.1 CASE 1 INPUT DATA 





In Figure 3.1, computed values for E[A,{t)] are plotted for 
several values of p. Simulation results are the actual mean values of aircraft 
availability at each t. Analytical solutions are computed at each time ¢ according to 


the formula 
E[A,(@)] = min {A,, K,-a.B,(@),..., K,-a.B,()}. (3.5) 


Formula (3.5) is correct to order a only; in practical terms it does not recognize the 
random variability of K;-N,(t). Note that, at all tf values examined, there is fairly good 


agreement between simulation and analytical solutions, even for small values of p. 
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Similar parameters 
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Figure 3.1 Case 1 Availability 


(2) Case 2. Table 3.2 shows that failure and repair rates for Case 
2 present much higher variation. As in Case 1, the sample system comprises 4.=50 
aircraft having each J=10 vital modules. All aircraft have exactly one module of each 


type installed. 
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Table 3.2 CASE 2 INPUT DATA 





As Figure 3.2 depicts, the choice of the parameter p now 
more strongly influences the accuracy of the analytical solutions. Note that p>30 is 
required to obtain a good approximation in this case. 

The importance of these initial results is that the user of the 
present analytical models should be aware that the selection of the appropriate 
tuning parameter is not case-independent. A reasonable implementation of these 
analytical models should consider a pre-verification of rate ranges before assigning 
the value of p. For technical reasons, however, it is desirable to make p as small as 


possible. 
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Figure 3.2 Case 2 Availability 


B. NORMALITY ANALYSIS 

An important result involving the diffusion approximation for large systems 
under heavy traffic conditions (1.2) is that N(t) in (2.2) is approximately multivariate 
normal with mean @f(t) and variance-covariance matrix aZ(t). Consequently, K-N(f) 
must be also normal with mean K - af(t) and the same variance. Here A(t) is the 


vector with components f(t), i=1,...,/. 
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Pilnick [Ref. 8, p.100-102] applied classical statistical analysis to the simulation 
data and demonstrated that the normality assumption appears to hold for particular 
transient times, as well as for the steady-state phase. In this section, alternative non- 
parametric methods are employed to verify the normality theory for the complete 
mission period. 

1. Methodology 

If the hypothesis of normality holds for A,{t), then the 5 and 95% 


percentiles of the distribution of A,{t) are approximately 


@o5 = EfAy(t)] + zo, x (Var[Ay()] 
@o5, = EfAy()] + 295 x ¥VarlAy@] , 


where z, is the a’ quantile of a standard normal distribution. 


(3.6) 


Values of E[A,{t)] and Var{A,{t)] can be computed by simulation and, 
importantly, also by an analytical-numerical method, and (3.6) can be used to 
calculate y,, and ~yo;. Now suppose samples of size n=1000 (i.e., 1000 replications 
in the simulation) are generated. If the 50" and 950" sample order statistics are 
stored for each time ¢, then, under the normal hypothesis, it is presumed that these 


values will approximate the theoretical numbers. 


2. Numerical Results 
The sample system represented by Case 2 above is used as a numerical 


example. Figure 3.3 displays the mean value of A,{t), the theoretical percentiles 
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computed by (3.6), and the corresponding sample order statistics. Note the jumps in 
the order statistics, denoted by A(S0) and A(950) in the plot. This is explained by the 
fact that they are integer values, representing actual aircraft availability at t for a 
particular replication. It is clear that the agreement is very satisfactory, confirming 
the usefulness of the normal model approximation. It thus becomes attractive to 
compute the probability that the number of available aircraft at time ¢ is less than 


or equal to x by use of a normal approximation. 


LAIN Policy 





Figure 3.3 Test for Normality 


C. CASE STUDY 
In Table 3.3 input data for a particular system are presented. This example is 


taken from Gaver, Isaacson and Pilnick [Ref. 9], who analyzed the effect of different 
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stockage patterns using the analytical formulation. Once more, failure and repair 
rates show large variation. The objective in this section is to reproduce the analysis 
of the relative performances of the adaptive (LAIN) and FCFS policies, using 
simulation and analytical techniques. 


Table 3.3. CASE STUDY INPUT DATA 
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Several modifications to the original set-up are introduced and the consequent 
effect on aircraft availability is analyzed. For all cases, simulation (S) and analytical 
(A) results are tabulated for both policies. These results are the mean values 
(E[AV]), standard deviation (SDEV), and the 5" and 95" percentiles (95% Cl) of 
the distribution of A,{t). Appendix A reproduces these results in graphic form; for 


each case that is analyzed in this section, plots are provided for comparative 
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performance for the two policies, as well as probabilistic limits for the actual 
availability using both simulation and analytical formulations. 

In this sample case there are initially A, = 50 aircraft. For the LAIN policy, the 
"anti-availability" parameter p is set to 30 in (3.3). Table 3.4 exhibits the resulting 


combat availability for each policy when no spares are provided, i.e., K;=50, 


1=1,2,...,10 (Case A). 


Table 3.4 COMPARATIVE PERFORMANCES (CASE A) 
FCFS-S | FCFS-A_ | LAIN-S | LAIN-A 
| | | 
t E[AV] SDEV | E[AV] SDEV | E[AV] SDEV | E[AV] SDEV 
| | | 
momeeceoy 9297 | 39.27 3.31 | 43.38 2.10 | 43.18 2.55 
PSD PS .2I6) | 32.76 3.78 | 40.22 2.06 | 40.31 2.51 
Bom 26866 3.63 | 28.81 3.89 | 38.01 2.13 | 38.13 2.40 
aome26-43 3.58 | 26.42 3.90 | 35.92 2.13 | 36.30 2.30 
Pomeeseio 63.65 || 24.95 3.90 | 34.36 2.10 | 34.71 2.21 
pomecamue 83543 | 24.08 3.89 | 33.03 1.96 | 33.32 2.14 
viemeocmay 3.60 | 23.54 3.89 | 31.81 1.85 | 32.10 2.07 
pemenoe ese O | 23.22 3.69 | 30.90 1.86 | 31.03 2.02 
emeeeeeo s.71 | 23.02 3.89 | 29.90 1.92 | 30.08 1.97 
HOOmMeeMSO0 4467 | 22.91 3.89 | 29.09 1.89 | 29.23 1.92 
| | | 
95% CI | 95% CI | 95% CI | 95% CI 
10 43.46 33.68] 44.72 33.82] 46.85 39.91] 47.39 38.97 
20 37.73 26.91] 38.99 26.53] 43.62 36.83] 44.45 36.18 
ome 34,65 22.67| 35.23  22.39/ 41.52 34.49] 42.09 34.17 
Aes Aue O15 132.6605 191991. 39.44 32.411 40.09 32.51 
50 31.21 19.16] 31.38 18.52] 37.83 30.89] 38.36 31.06 
60 29.72 18.39| 30.50 17.65] 36.26 29.81] 36.85 29.79 
70 29.40 17.54] 29.96 17.13] 34.86 28.77| 35.52 28.68 
80 29.03 16.82] 29.63 16.81] 33.96 27.83| 34.35 27.70 
90 29.12 16.87] 29.44 16.61] 33.06 26.73] 33.32 26.83 
100 28.55 16.45| 29.32 16.49] 32.04 26.13] 32.40 26.06 
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As it was expected, the LAIN priority policy achieves larger average values for 
combat availability at all times. Note, also, that the adaptive scheme induces lower 
variances than does the simple FCFS policy. It should be noted especially that the 
diffusion approximation results match those from the simulation quite accurately. At 
least for this example the standard deviations derived using analytical method (A) are 
slightly larger than those from simulation (S). If this slightly conservative bias prevails 
it is a desirable situation. 

In Case B one attempts to increase combat availability by providing two spares 


for each module. Table 3.5 displays the results for Case B. 
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Table 3.5 COMPARATIVE PERFORMANCES (CASE B; 2 SPARES) 
ECES=S | FCFS-A | LAIN<-S | LAIN-A 
| | | 
1 E{AV] SDEV | E[AV] SDEV | E[AV] SDEV | E[AV] SDEV 
| | | 
CMO So. lye| 40853053 -398 | 45203, 1.96] 44.65, 2.65 
meme 5S. 67 imo omozm 3.86 | 41-48 2.02 | 41.51 2.60 
Sepeetiowso 3-5) | 292805 3.97 | 38.96 2.23 | 39.17 2.49 
Come cigeaoems.7/ | 2671) 3.97 | 36.75 2.18 | 37.22 2.38 
pm accor 3.84) 25.14 3.95 | 35.174 2.17 | 35.53 2.28 
SOmezames 3-96 | 24419 3.94 | 33:59 2.06 | 34.06 2.20 
Wome comes 4.08 | 23462 3.94 If 32:41 2.00 | 32.76 2.12 
Pmeeeore oe 3. | C2 Oe hCU8Ul ISP 6US12306CU«dLCWC OS} 3166206 62.06 
Pumeevomeces S.06 | 23205 3.93) | 30.37 1.89 | 30.62 2.01 
iemecomese  SeGD | 22.92 3.93 | 29.55 1.92 | 29.73 1.96 
| | | 
95% CI | 95% CI | 95% Cl | yo5 Cl 
| | | 
10 45.11 34.66] 46.11 34.94| 48.26 41.80] 49.02 40.29 
Pes InoOOn 26.59 | 39.89 27.-16| 44.81 38.16] 45.81 37.22 
iOmesemtom 22.0600) 35.82 @22.73| 42.64 35.29| 43.29 35.06 
come sem675 20.247 33.26— 20.15| 40.35 33.16| 41.14 33.30 
pemees teoeme be.w6o6) 31266, 18.62) 38°73 31.56|/) 39.29 31.77 
TOMO lo.2zo) 30269" 17.70| 362.99 30220] 37.68 30.43 
pees Oro Lo.o4) 30.11) 72.12) 35271 29710) 36.27 29.26 
Wome or 16.671) 29275 W6e.78]| 34241 28.05] 35.03 28.22 
See oodee 16.92) 29.54 6.56] 33.49 27.24] 33.93 27.31 
MOO Zoi? | 16590) 29.41 16.43) 32.71 26.39| 32.96 26.49 


The numbers show a small improvement in the average availability, specially 
in early periods; however, there is very small improvement at later times. 
The effect of providing ten spare parts for each module (Case C) is examined 


next; see Table 3.6. 
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Table 3.6 COMPARATIVE PERFORMANCES (CASE C; 10 SPARES) 
FCFS-S_ | FCFS-A_ | LAIN-S | LAIN-A 
| | | 
t E[AV] SDEV | E[AV] SDEV | E[AV] SDEV | E[AV] SDEV 
| | | 
10 45.56 3.23 | 46.40 3.83 | 49.53 1.80 | 50.00 3.42 
20 36.84 3.97 | 37.17 4.19 | 46.21 2.69 | 46.61 3.26 
30 31.17 4.01 | 31.46 4.24 | 43927 2.809) 43e@) 3.08 
40 27.93 3.93 | 28.01 4.20 | 40.62 2.65 | 40.83 2.87 
50 26.03 3.90 | 25.93 4.15 | 318.47092.54 | 3880 =2mec 
60 24.74 4.07 | 24.66 4.12 | 36.65 2.48 | 37.03 2.51 
70 23.67 3.97 | 23.89 4.10 | 35.16. 2.36 [e35ee5ue2 0 
80 23.42 4.02 | 23.43 4.09 | 33878 2.30 |B S4mocumc em 
90 23.06 3.77 | 23.15 4.08 | 3262. 2.220) 32eeo eee eee 
100 23.01 3.83 | 22.98 4.08 | 31.68 2.108] 315745026 
| | | 
95% CI | 95% CI | 95% CI | 95% CI 
| | | 
10 50.89 40.24] 52.91 40.28] 51.50 45.55] 55.65 44.35 
20 43.39 30.28] 44.09 30.26| 50.66 41.77] 51.99 41.22 
30 37.67 24.46] 38.46 24.46| 47.77 38.73| 48.38 .21 
40 34.41 21.45] 34.93 21.09] 44.96 36.23] 45.56 36.09 
50 32.57 20.02] 32777 19.08| 42256) 34388 4359 same 
60’ 31.47 18.34] 31.45 17.87] 40.70) 3255018 41.07 oomee 
70 30.04 17.30] 30.65 17.13] 38.98 31.20] 39.41 31.49 
80 29.96 16.75| 30.17 16.68] 37.65 29.91] 37.86 30.25 
90 29.28 16.84] 29.88 16.41] 36.28 .95| 36.50 29.17 
100 29.34 16.68] 29.71 16.25| 35.14 35.30 28.19 


28.22] 


A significantly greater effect on aircraft availability occurs in Case C. The 
improvement in availability is specially noticeable in early times for both policies, 
with the effect decreasing considerably at the end of the mission period. Note, also, 
that the effect is more important in the LAIN availability. In the LAIN case, there 
is a much larger relative error between the diffusion approximation and simulation 
with respect to standard deviations for early times; again the standard deviation 


obtained analytically (LAIN-A) noticeably exceeds the simulation value (LAIN-S). 
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It is interesting to observe what happens when ten spares are allocated only to 


those five modules with higher failure rates (Case D); see Table 3.7. 


Table 3.7 COMPARATIVE PERFORMANCES (CASE D) 
FCFS-S_ | FCFS-A_ | LAIN-S- | LAIN-A 
| | | 
t E[AV]) SDEV | E[AV] SDEV | E[AV] SDEV | E[AV]  SDEV 
| | | 
momer46505002.44 | 46.99  3.98.| 46.17. 1.83 | 46.11. 2.90 
pemees Gees. 78° | 37.43) 4.21° | 42.72) 2.07 | 42.59 2.74 
pOmeenens 24-27 sl 31,.6204.24)) 39.80 2.19.1 39.97. 2.63 
Orme 05S. 90 nile es om 4.20m| 37.55 2.23 | 37.95) 2.50 
HOM 26002 3.67"9| 25.98 4.15 | 35.62 2.30 | 36.22 2.36 
pOmeedeyameed.O8 | 24.699 4.12 | 34514 2.077] 34.64 2.24 
Ome erciee2 95 | 23.90" 4.10,| 32.91 1.999] 33.24 2.17 
perc aes.99 |) 23.5448 4.08 | 31.87 1.937] 32.01 2.09 
Bommeeeetoe s.c4 | 23.16 4.08 | 30.79 1.91 | 30.94 2.03 
MoOmese com a.Se | 22.98 4.06 | 29484 2.00 | 29.99 1.98 
| | | 
95% CI | 95% CI | 95% CI | 95% CI 
| | | 
10 48.97 40.92] 53.56 40.43] 49.18 43.16] 50.89 41.33 
20 42.98 30.50] 44.37 30.49] 46.13 39.31] 47.11 38.07 
piomee 6 24.23) 38.61 24.62). 43.42. 36.18| 44.31 35.62 
Aomeo4e4s 21.62] 35.03 21.18) 41.24 33.87| 42.08 33.82 
pOmeal tT 20.071) 32.83 192513] 39.41 31.84/940.11 32.32 
omen 15.12) 31.466 17.90| 37.55 930.73) 38.34 30.95 
Fipme sons oe 17.30| 30.66" 17.14] 36.19  29.62| 36.80 291.67 
80 30.13 16.96] 30.18 16.70] 35.06 28.69] 35.46 28.56 
Bon) 296530) 17.16| 29.89 16.42| 33195)° 27.64] 34.29 27.59 
mMOON29= 55 16.90| 29.71 16.26| 33.14 26.53] 33.26 .73 


In this case, the FCFS policy performs slightly better than the priority scheme 
in the beginning of the mission. The implication may be that an appropriate 
weighting, taking into account high failure rates, may prove to be more efficient than 
the exclusive concern with the current least available module at each time. The 


LAIN policy is again more efficient for later periods. 


29 


Now suppose that all modules cost the same (not very realistic, of course), and 


instead of allocating ten spares to those modules with higher failure rates, as in Case 


D, five spares are introduced across the board (Case E); see Table 3.8. 


Table 3.8 COMPARATIVE PERFORMANCES (CASE E; 5 SPARES) 
FCFS-S | FCFS-A_ | LAIN-S_ | LAIN-A 
| | | 
t E[AV] SDEV | E[AV] SDEV | E[AV] SDEV | E[AV] SDEV 
| | | 
10 42.00 3.22 | 42.59 3.51 | 46.95 2.70 | 46.95 2.83 
20 34/14 3.70 | 34.75 3.97 | 4363602. Som) 4cmcecmo wee 
30 29.63 3.85 | 30.00 4.06 | 40.47 2.40 | 40.60 2.66 
40 27.03 3.80 | 27-12 4.04) || 38320 92. 29m) somo em 
50 25.29 3.81 | 25.38 4.08 | 36.82 82.206) 36671 oeeo 
60 24.27 3.86 | 24.33 4.01 | 34992 $2.20) somo so 
70 23.68 3.87 | 23.69 4.00%] 33944 72.106) couvc cee 
80 23.43 4.03 | 23.32 3.99 | 32030 2030) g2y5 ieee 
90 23.27 3.79 | 23.09 3.99 | 31.19 1598 | 3242) 2eos 
100 22.99 3.66 | 22.96 3.99 | 30.22 202 1) 0m. Wee 
| | | 
95% CI | 95% CI | 95% CI | 95% CI 
| | | 
10 47.63 36.71] 48.39 36.79| 50.36 43.54] 51.63 42.28 
20 41.20 28.07] 41.32 28.19] 47.38 39-26] 47-92eocme 
30 35.87 23.38| 36.71 23.29] 44.43 36.51] 45.10 36.30 
AQ 33.11" 20.96] 33382 20.43] 42¥00" 34.49)) 4257 46a 
50 31.50 19.07] 32.04 18.741) 39.95) .02960| 4em77meoomad 
60 30.64 17.89] 30.95 17.72] 38.46 31.18] 38.95 31.37 
70 30.07 17230| 30880 17.120] 36.90 29ecc)\" 37 et CcOmnO 
80 30.09 16.80] 29.91 16.72] 35.65 28.95] 36.06 28.98 
90 29.53 17.02] 29.67 16.49| 34.45 27.93] 34.86 28.00 
100 29.03 16.95] 29.53 16.35| 33.55 26.89] 33.81 27.13 


It is illustrative to observe that this new, and almost certainly less expensive, 
stock plan influences each scheduling policy in a different way. In the FCFS case the 


combat availability declines throughout the mission when compared to Case D, 
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especially at earlier times; on the other hand, and somewhat surprisingly, LAIN 
availability is slightly enhanced at all times. 

Table 3.9 shows the result of submitting the system in Case E to an arbitrary 
weighting consisting of w,=v,/), i.e., a high weight for modules with high repair rate 


and low failure rate, using the LAIN policy. 


Table 3.9 THE EFFECT OF WEIGHTS IN CASE E 


WEIGHTS NOT USED WEIGHTS USED 


| 46.9 
43,3 


Shi 
31.4 
30.4 
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Quite surprisingly, the effect of such weighting procedure is insignificant (in 
fact, combat availability decreases slightly). It is necessary to understand that the 
numbers in Table 3.9 were generated by setting w,=7 and w,=v,/,, 1=1,2,...,10, in 


(3.2). The result is obviously the same if the w; are set to y,/; instead of v,/),. It is 
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not clear if this fact implies that this arbitrary weighting is ineffective, or if the model 
does not truly represent such strategy. 

Finally, simulations using various weighting factors were run for the LAIN case, 
including the procedure in which the module type selected for repair next is that with 
the least value of the product of the actual availability by the ratio failure/repair 
rates. All the results showed that the combat availability is reduced, when compared 
to the unweighted LAIN scheme. 

In our last sample system (Case F) the parameters are modified in such a way 
that the heavy traffic condition (1.2) for the diffusion approximation is barely 


satisfied; see Table 3.10. 


Table 3.10 CASE F; INPUT DATA 


MODULE 
1 
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Table 3.11 shows the computed expected values and standard deviations for A,{t) for 


If (1.2) is applied to the data we have 


10 
7 
SUE ee = aan . 
i=] v; 


both the simulation (S) and the diffusion approximation (A). 


Table 3.11 COMPARATIVE PERFORMANCES (CASE F) 
FCFS-S FCFS-A LAIN-S LAIN-A 

t E[AV] SDEV  E[AV}) SDEV  E[AV] SDEV E[AV] SDEV 
Mame4gueo 2.625 44.28 2.73 46.58 1.46 46.22 2.80 
Pome cio 3.30 41.93 3.34) 45.17) 1.68 44.97" 2.15 
BlOmes9650)3-64 40.21 3.53 944906 1.77 43.91 2.10 
Moe See60 3.80 38.94 3.59 43.09 1.87 43.00 1.85 
SS Oe eo Nl i ADO G eS 6 MAD. 21 2 327 
60 37.51 3.65 37.96 3.66 41.54 1.96 41.52 2.20 
moms?) 3.65.037.73.3.77. .40.94 1.95 40.91 2.52 
80 37.07 3.75 37.44 3.92 40.35 1.98 40.39 1.99 
gmoye1O ) 3.64 37.25 3.87 39.82 1.98 39.93 2.30 
OMe eed ONO 713,44 0 9A DD 10 1 39. 53. 25 


It can be noticed that the approximation provided by the diffusion model is still 
very useful, especially for mean values; the (now larger) error in standard deviation 
is, again, biased upward. A case where the heavy traffic condition is violated is 


investigated in Appendix D. 
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The effect of supplying ten spares for each module (Case G) is examined in 
Table 3.12. Once more, the effect is more important in the LAIN availability, 


especially for t>60. 


Table 3.12 COMPARATIVE PERFORMANCES (CASE G; 10 SPARES) 

FCFS-S FCFS-A LAIN-S LAIN-A 
t E[AV] SDEV  E[AV] SDEV  E[AV] SDEV’ E[AV]  SDEV 
10.49.50 1.27 ...50.00°92.43... 49.62 O7CmD Oo OO elec 
20 46.44 3.36 °&-47.00 4.16 49-54 0.86 "Senco iD 
30m 43.64 3.84 44.10 004.274 OnGie 1. 22 ono ome eo 
40 41.04 3.98 41.83 4.71 48.32 1.60 49.36 1.83 
50 39.49 4.11 40.17 4.89 47839. 1.5930 47e76umoeeo 
60 38°57 4.10 39.07 4.49% 46794499 2.09) 4637 ue2ees 
70 37.94 4.22 38.52 4.47 45.38" 2.159) 4550c oa 
80 37.51 4.05 38.03 4.239 44.417 2.2065 44 ome oes 
90 37.40 4.06 3757379 5.05 43.435 2.295 43 eicmeomee 
100 37.10 4.02 37.44 94.6195 420670. 26m 4256c meee 
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IV. FURTHER DEVELOPMENTS 


A. NON-CANNIBALIZATION 

In all models considered so far it is assumed that working components 
belonging to inoperative aircraft (resulting from previous failures) can be used to 
substitute for modules which have just failed in another (otherwise operative) 
aircraft. This circumstance is what is known as "cannibalization", and it is usually 
implemented in real situations when it is necessary to maintain a high level of 
combat readiness through increased availability of vital assets for a limited period of 
time. As an example, the overall objective of the Chief of Naval Operations is to 
obtain at least seventy-two percent of fully-mission-capable aircraft in a squadron 
(Ref. 1:p.4]. 

In this section a modification of the original "analytic" model is developed for 
a situation where cannibalization is not employed. It is assumed that the time horizon 
of interest is relatively small, such that the adoption of a cannibalization policy does 
not affect significantly components’ failure rates’. 

It is further assumed that an aircraft will return to its base almost as soon as 


the first vital component fails, and that no other component in this aircraft will suffer 


‘Experience shows that this kind of maintenance policy tends to decrease 
equipment life time if utilized for long periods. The effect is related to an increase of 
the individual failure rates, A; This fact is corroborated by this author’s own 
experience as head of the electronic equipments maintenance division in a Brazilian 
Navy FFG (Guided Missile Frigate). 
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failure until it is back in operation. This is, admittedly, a somewhat unrealistic 
assumption, but the lower bound it represents may be useful to demonstrate the 
dramatic role cannibalization plays in terms of increasing availability. The model 
may, also, serve as a quick tool for assessing the effect of this kind of maintenance 
policy. A model that supposes that degraded aircraft proceed with their mission can 
also be constructed, but this is not done at this time. 

Let A,,.(t) denote the total number of operational aircraft at each time ¢, for the 
non-cannibalization case. At t=0, let A,.(0)=A,, the number of initially deployed 
aircraft. 

As is the case now, a failed component renders an aircraft inoperative until it 
is repaired, assuming there is no spare part in stock. Otherwise, the failed component 
is Substituted, and will join the appropriate queue. Let S; = K; - A, be the number 
of spare parts available for module 7. 


A,,-(t) may be represented by 


A, (ft) = A, - : T(t)x[ Nd - S,] (4.1) 


where /, K;, and N,(t) are as defined in chapter I, and /,(t) is the indicator 


variable defined by 


He INO) SS. (4.2) 
MO = 10 if Nit) <S, 
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1. A Numerical Example 

Once again, a Monte Carlo simulation model is used to validate analytical 
results. The input data for the example used in this investigation are those of Table 
3.2 in Chapter III. 

In the simulation model, A,,(f) is updated only when a failure occurs or 
when a repair is completed, after considering inventory levels. When a failure occurs, 
the aircraft is brought back to the operational status if a spare exists, and the value 
of A,(f) remains unchanged. Otherwise, the aircraft must wait until repair 
completion, and A,,(t) decreases by one. When a module is repaired, 4,,-(f) increases 
by one only if the component was causing an aircraft to remain inoperative. 

With respect to the equations for Bt) in (2.4), the only necessary 
modification is the use of A,,(t) (4.1) scaled by a in the place of a,{t). For the non- 
cannibalization case, only the LAIN policy is examined, and the corresponding results 
are compared to the usual (perfect cannibalization) plan, for both FCFS and LAIN. 

The numerical example attempts to exploit a situation where the modules 
show large diversity in terms of failure and repair rates, and p=30 is used in (3.4). 

Figure 4.1 displays a fairly satisfactory agreement between simulation and 
"analytical" results. Numerical values for the expected value of the aircraft availability 


as a function of t, te[0,100] are shown for both cases. 
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Figure 4.1 Non-Cannibalization LAIN Policy 


Figure 4.2 exhibits the comparison among the various policies, showing the 


evident decrease in availability when cannibalization is not employed. 
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Comparative Performance 


E[Avi(t)] 





Figure 4.2 Non-Cannibalization Performance 


Table 4.1 summarizes these results for t= 10,20,...,100. The numbers are 


rounded to the closest integer. 
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Table 4.1 NON-CANNIBALIZATION PERFORMANCE 





B. EXTERNAL REPAIR (BEYOND THE CAPABILITY OF (LOCAL) 
MAINTENANCE - BCM) 

A more realistic characterization of the Air Unit -Repairman system must 
consider the case in which failed modules require a degree of support that is beyond 
the capability of local maintenance. The implication is that a certain portion of failed 
modules must be sent to a higher maintenance level, returning to the local level after 


an uncertain time. See Gaver, Isaacson and Pilnick [Ref. 9]. 


1. Pre-Local-Repair BCM 
Suppose that every time a module of type j fails there is a probability n; 
that the local repair shop will not be able to service it. This implies that this module 


has to be sent to another repair echelon, external to the deployed group. Denote 
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these modules as BCM modules. Further assume that all BCM modules, after a 
random period of time, will return to the system. 

The amount of time a BCM module of type j stays outside the system (to 
be repaired / substituted) is assumed to be an exponential random variable with 
mean 1/6,, i.e. BCM —, type j return to the local level at a rate equal to 6,. 
Let K,(t) be the time-dependent number of modules of type j available at the local 
level at time ¢, with K,(0)=K,, K;-K,(t) BCM modules at time ¢, j=1,2,...,/. 

For all j=1,..,/, the probability P;(N(¢)) = P,(¢) that a module of type j fails 
in the interval (t,+dt) and is repairable at the local repair facility (non-BCM), 


conditional on the present state N(f) is given by 
P(t)dt = (1 - n,)4,[K\@) -Nj@)] dt + o(at) (4.3) 
Now define a(t) ~ K;(t)/a, and substitute a(t) for a, in (2.5) for the 
definition of a’{t). 
The following systems of differential equations may be written for the 


approximate scaled mean of N(t), B(t), and for the scaled number of modules at the 


local level a(t): 


dB (t) 








dt ACL - 1; )a'yO) - HP) (4.4) 
da (t) : | 
dt 6 (a; a(¢)) - A;nja’ At) 
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a. A Numencal Example 

Consider the sample system depicted in Chapter III-C (Table 3.3), and 
suppose that all BCM modules return to the local level at a rate 6;= 1/30, i=1,2,...,10. 
Further assume that the proportion of BCM modules of type j, 7;, is equal to 0.20 for 
all i. Figure 4.3 displays the mean aircraft availability for te[0,100], calculated 
according to (4.4) and (4.5) for both the FCFS and LAIN policies. The plot includes 
simulation results, and it is clear that the accuracy of the analytical model is very 
acceptable for this case. 

The correctness of the model was tested for different values of 7;. 
However, numerical and analytical difficulties are observed for some particular 
parameter values. A possible explanation may be that the decrease in local demand 
weakens the heavy traffic condition necessary for the present formulation. In 


Appendix B this subject is further explored. 


42 


E[Avi(t)] 





Figure 4.3. Pre-Local-Repair BCM 


2.  Post-Local-Repair BCM 
Gaver, Isaacson and Pilnick [Ref. 9] conceived an interesting modification 
to the pre-local-repair BCM case. They argue that immediate consignment to 
external repair of a failed module is, to say the least, an optimistic assumption. If 
every failed module must be first submitted to the local test and repair cycle, and 
only at termination it is either completely repaired or a decision for external repair 


is made, they model the modified system using 
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aB (t) 
— = daft) - 1p) (4.6) 
da. 

i = §(@,-a,(t)) - 1jn;P,(t) Goll) 





a. A Numerical Example 


Once more consider the system described in Table 3.3, and suppose 
that 6; = 30 for all i. Figure 4.4 and 4.5 show a useful accurate agreement between 


simulation and analytical results for different values of n,. 
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Figure 4.4 Post-Local-Repair BCM (n = 0.20) 
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For the post-local-repair case, numerical and/or analytical 
complications were not observed when the model was tested under a considerable 


number of parameter combinations. 


eta = 0.40 


E[Av(t)] 





Figure 4.5 Post-Local-Repair BCM (n = 0.40) 
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V. OPTIMIZATION MODELS 


In this chapter two simple optimization models are described for the perfect 
cannibalization, non-BCM case. The models attempt to represent different decision 
problems, and can be readily adapted for different situations, such as BCM. One of 
the basic assumptions here is that an adaptive scheduling policy (LAIN) is used. 

Every stochastic optimization problem must reflect an attitude toward risk; it 
is therefore often possible to pass from one type of formulation to an equivalent one. 
The models described here are designed for planning purposes, in the sense that the 
decision does not depend in any way on future observations of the random vector 
N(t). Such models are usually called anticipative models; see, for example, R. J.-B. 


Wets [Ref. 14]. 


A. DEFINITIONS 
In order to describe the optimization models for the combat availability 
problem, a larger set of variables must be defined. The complete set of variables is 


defined below: 


Let 
I = number of distinct mission-essential component types in each a/c; 
1S, = initial provision of components of type /, including those installed 


in the aircraft ; 
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C; = cost of one unit of component type ; 
A. = number of aircraft initially deployed ; 
S; = spare parts provided for component of type i, i.e., $;=K;-A, ; 


time horizon ; 


i 
P.»in = Yinimum percentage of the original number of a/c required at time 


s(t) = argmin {K,(t) - N,(t)} ; 

Aft) = min {A,, K-N,(t)} ; 

Ahgee = 

N,(t) = number of components of type 7 being repaired or waiting in 


queue at time ¢ ; 


B,(t) = approximation of the scaled mean of N,(t) ; 
2(t) = scaled variance-covariance matrix of N(t) , and 
a = A, + 2; K; . 


As was shown in chapter III, N,(t) is approximately normally distributed, 1.e., 


N@) = N(a.p(),a.2)) . (5.1) 


where N(f) is the vector with components N,(t), and B(t) has components £(f), 


een ae 
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B. RELIABILITY MODEL 

Suppose that, during the planning phase of a combat mission that is expected 
to last T units of time, it has been established that a minimum percentage P,,,,, of the 
initial number of fully mission capable aircraft is required with high probability, r. 
Then, assuming that resupply is not possible during the mission, it seems reasonable 
to consider only the combat availability at time 7. Further imagine that a general 
stockage policy for many missions is under study, so that costs must be minimized for 
this particular mission. 


A general mathematical description of this problem is 


min > ¢,S; (7) 
SP lA) 24 he ee el (ii) () 
5 senna ea) 


We can interpret (5.2) as a wish to minimize spares allocation cost subject to 
a minimum system reliability, where the system is now viewed as the whole Air Unit. 
In order to correctly employ the approximation in (5.1), it is necessary to 
examine the probabilistic (or chance) constraint (i) in (5.2) with respect to the 


corresponding expression for each item. Consider 


AT) = min{A,,K, -N(T)} 


min {A,,(4,+S,) -N{T)} 


min{4¢, (Ag+ Syn) - Nyy (T)) 
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Now, under heavy traffic and for large enough 7, so that combat availability is 


determined by the module with the least availability, i.e., 


AT) = Ac + S5(1) i. Nn) ) 
the event A,(T) > Av,,,, can be rewritten as 
(A+ Sr) - NenP) > Ay | 
2 NnF) < (4 -+ s(T)) Aviad 


and, using (5.1), we have approximately 


(A, + Ss) ~ Ay, ~ FB sen) 


P{N < (A.+8 -A a 
Bsn) $ Ae * Seen) Arua] 4.0 .n(T) 





where 0,(7) is the diagonal element of the matrix 2(T) corresponding to the 
module s(7). 


A deterministic approximation for the model in (5.2) is, then 


min JY c;,S, 


St. Soop 2 2,40 .)(T) + Ay + @B p(T) - A, (5.3) 


S.>0 : 


l 


where z, is the ”” quantile of the standard normal distribution. 
Note that (5.3) is a continuous approximation of the actual discrete problem 
(5.2), and that, even though the constraint in (5.3) is described with a single 


inequality, in practice its computation at each stage of the optimization procedure 
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involves the solution of the system of differential equations used to define A(t) and 


MG) 


C. EXPECTED VALUE MODEL 

It is common practice in stochastic programming models to make use of 
expected values when the coefficients of the decision variables are random. In the 
following simple model we seek to achieve the maximum expected value of combat 


availability at T, subject to a budget constraint: 


max E[A,{T)] 
I 

co DL C DD (5.4) 
i=1 
S COM ese 


where B is the budget. 


The approximation to the objective function in (5.4) 


E[A,(7)] = mint4,,(A,+S,) - a8 ,(7),...,(A.+S,) -aB (1) 


may,once again, be simplified to 


A. + S51) 7 aB 7) 


under the same conditions described in the last section. 
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The continuous approximation to the problem (5.4) is, then 
max (A, +S.) - 4B p(T) 
I 
SSS GS, 7a (5.5) 
i=1 


Saa0 


i 


ol 


VI. CONCLUSIONS AND RECOMMENDATIONS 


In this thesis several mathematical models for evaluating scheduling and spares 
stockage rules in a transient dynamic combat environment are reviewed, developed 
and analyzed for a situation in which a large population of modules under heavy 
traffic conditions drives combat availability. Even though this paper specifically 
considered an Air Unit detachment problem, it is clear that these models may be 
utilized in numerous environments. 

The adaptive priority policy (LAIN) was shown to considerably improve aircraft 
availability when compared to a standard FCFS scheme. The modeling technique 
using a diffusion approximation provides speedy solutions for a wide range of 
problems, and it provides the means for ready comparison of alternative scheduling 
policies, reflecting diverse organizational maintenance disciplines. 

Initial results suggest that a pure priority scheme, based on the least available 
item at each time, overperforms arbitrary weighting procedures. This subject must 
be further explored. 

The model for situations where cannibalization is not employed may help a 
decision maker, through a quick demonstration of the dramatic effect of this policy 
in terms of combat availability. 

The situation in which failed modules require a degree of support that is 


beyond the capability of the local echelon is readily accounted for with small 
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modifications to the original model. In this thesis we describe these modifications for 
the mean availability only. An area for further research is the development of 
equations for the variance of these BCM systems. Another point of interest is the 
development of models that take into account some mixture of pre- and post-local- 
repair allocation to distant — as well as possible mistakes in these consignments. 

The optimization models described here are an attempt to apply the analytical 
models as a framework for choosing spare modules allocations in various 
environments. The actual solution of these models is certainly a challenging 


continuation of the present work. 
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APPENDIX A. GRAPHICS FOR THE CASE STUDY 


In this appendix, a series of plots are used to demonstrate graphically the 
precision that is achieved by the diffusion approximation models. Cases A through 
D in Chapter III-C are reproduced here. For each case, the following plots are 
provided: 

1. FCFS vs. LAIN - the expected combat availability for both policies is 
graphed for te[0,100]. Simulation results are also plotted, for comparison; 

2. FCFS policy - lower and upper probability levels are shown, together with 
the mean availability. These levels are computed as the expected value of combat 
availability (E[A,y(t)]) + the standard deviation (s) times the 95" percentile of the 
standard normal distribution (1.65). Simulation values are also plotted; 


3. LAIN policy - as in item 2 above. 
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FC-FS vs. LAIN 
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Figure Al. Case A; Both Policies 
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FC-FS (Case A) 





Figure A2. Case A; FCFS 
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Figure A3. Case A; LAIN 
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Figure A4. Case B; Both Policies 
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FC-FS (Case B) 
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Figure A5. Case B; FCFS 
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LAIN (Case B) 
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Figure A6. Case B; LAIN 
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Figure A7. Case C; Both Policies 
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FC-FS (Case C) 
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Figure A8. Case C; FCFS 
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Figure A9. Case C; LAIN 
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FC-FS vs. LAIN 
Case D 
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Figure Al0. Case D; Both Policies 
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FC-FS (Case D) 
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Figure All. Case D; FCFS 
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LAIN (Case D) 
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Figure Al2. Case D; LAIN 
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APPENDIX B. PRE-LOCAL-REPAIR BCM 


In this appendix, some of the numerical/analytical difficulties that occur in the 
pre-local-repair BCM model in (4.4) and (4.5) are summarized graphically. The 
model was extensively tested by means of a very accurate ODE solver developed by 
W. B. Gragg (Naval Postgraduate School) for the MATLAB-PC package [Ref. 13]. 

The model was tested for several combinations of values for the parameters 7 
and 6. Figure 4.3 in Chapter IV displayed very accurate results when n;=0.20 and 6;7 
= 30, 7=1,2,...,10. For some parameter values, however, even stiff ODE solvers fail 
to achieve convergence at certain time intervals. 

Figure B1 shows the result of applying our model to the sample system 
characterized by Table 3.3 in Chapter III, with n;=0.30 and 0.40, respectively. It is 
clear that large discrepancies exist, even at early times, for both FCFS and LAIN 
policies. A critical singularity at t=42 forces the approximation to diverge strongly for 
the FCFS availability when n=0.40. 

We may conclude that a modification to the original model (4.4) and (4.5) is 
necessary to account for these difficulties. We do not attempt to investigate these 


modifications in this study. 
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Pre - local - repair BOM 
eta = 0.30 





Pre - local - repair BCM 
eta = 0.40 





Figure Bl. Pre-Local-Repair BCM Divergence 
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APPENDIX C. DIFFUSION APPROXIMATION ROUTINE 


SUBROUTINE DIFFUS (NTIME,MAXTIM,AVERAG,STDDEV, 
& DMEAN,DSTD,ITYPE) 

INTEGER QUANT 

PARAMETER (QUANT = 10) 

INTEGER NTIME,MAXTIM,ITYPE 

REAL AVERAG(NTIME),STDDEV(NTIME) 

REAL DMEAN(NTIME,QUANT),DSTD(NTIME,QUANT) 


, MODEL VARIABLES 

INTEGER ONREP,DOWN(QUANT),K(QUANT) 
INTEGER NUMMIN,AGP,AV,TOTQU,AVMIN 

REAL LB(QUANT),NU(QUANT),W(QUANT) 
COMMON/MODEL1/ONREP,DOWN,K,NUMMIN,AC, 
& P,AV,TOTQU,AVMIN 
COMMON/MODEL2/LB,NU,W 


INTEGER INOP(QUANT) 
COMMON/TIME0/INOP 


er ee POLICY VARIABLE 
INTEGER POLTYP 
COMMON/POLICY/POLTYP 


i i DIFFUSION VARIABLES 
REAL MU(QUANT),ALPH(QUANT), Y0(QUANT),ACNORM,X0 
COMMON/APPROX/MU,ALPH, Y0,ACNORM,X0 


ee Te LOCAL VARIABLES 
INTEGER I,J, TSTEP 
REAL A 


KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKEKKEKKKKKKKKKEKKKKKKKEKKKEKKEKEKEKKEKEKE ES 


MAIN PROGRAM 7 


KKKKKKKKKKKKKEKEKEKEKEKKKKKKEKKKKKKKKKKKKEKKKKKKKKEKEKEKKEKHKEKKEKEKKKEEKEKEKS 


Pe eee SCALING OF INPUT PARAMETERS 
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A = REAL(AC) 
DO 10 I=1,QUANT 
A =A + REAL{(K(I)) 
10 CONTINUE 
ACNORM = REAL(AC) / A 
POLTYP = ITYPE 
DO 20 I=1,QUANT 
ALPH(I) = REAL(K(I)) / A 
MU(I) = NU(I) /A 
YO(I) = REAL(INOP(I)) / A 
20 CONTINUE 
> | ie... oe COMPUTE TIME STEP 
TSTEP = MAXTIM/NTIME 


SoeeeeM ee se ssnateuee CALL ODE SOLVER AND RETURN THE VALUES 
OF THE MEAN AND STANDARD DEVIATION FOR 
COMPONENT WITH LEAST AVAILABILITY AT 
EACH TIME, FOR TYPE "ITYPE" 

CALL DFMEAN(NTIME,MAXTIM,AVERAG,STDDEV,A,DMEAN,DSTD) 

RETURN 

END 


* * *% 


KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK KKK KKK KKKKK KKK KKK KKK KK KK KKK K K 


SUBROUTINE DFMEAN(NTIME,MAXTIM,AVERAG, 
& STDDEV,A,DMEAN,DSTD) 


KKKKKKKKKKKKKKKKKKKKKK KKK KKK KKK KK KKK KK KK KK KKK KKK KKK KKKKKEKKKKKKKKKKK EK 


et LOCAL ENVIRONMENT 
INTEGER QUANT,MXPARM,MXSTEP,METHOD 
REAL TOL,HINIT 
PARAMETER ( QUANT = 10 ) 
PARAMETER ( MXPARM = 50_ ) 
PARAMETER ( MXSTEP = 3000 ) 
PARAMETER ( METHOD = 1 __ ) 
PARAMETER (TOL = 1.0E-5) 
PARAMETER (HINIT = 1.0E-3) 


af ee INPUT/OUTPUT PARAMETERS 
INTEGER NTIME,MAXTIM 
REAL AVERAG(NTIME),STDDEV(NTIME),A 
REAL DMEAN(NTIME,QUANT),DSTD(NTIME,QUANT) 


ry MODEL VARIABLES 
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*%* *% * *%* * 


INTEGER ONREP,DOWN(QUANT),K(QUANT) 
INTEGER NUMMIN,AC,P,AV,TOTQU,AVMIN 
REAL LB(QUANT),NU(QUANT),W(QUANT) 
COMMON/MODEL1/ONREP,DOWN,K,NUMMIN,AC 

& P,AV,TOTQU,AVMIN 
COMMON /MODEL2/LB,NU,W 


i.e DIFFUSION VARIABLES 
REAL MU(QUANT),ALPH(QUANT), Y0(QUANT),ACNORM,X0 
COMMON/APPROX/MU,ALPH, Y0,ACNORM,X0 


eee, LOCAL VARIABLES 

INTEGER IDO, ISTEP, I,J, INTIME,IMIN 

REAL AMAT(1,1),PARAM(MXPARM) 

REAL FCN,FCNJ 

REAL X,XEND,OK(QUANT),AVNOW, Y(QUANT + QUANT* QUANT) 


EXTERNAL FCN,FCNJ,.IVPAG,SSET,ISMIN 


y(1), y(2), .... Y(quant) represent actual values for the scaled means. 
Y(quant+ 1), y(quant+2),..., Y(quant+quant*quant) are a vector representation 
of the variance-covariance matrix in column form. 


Bere eee Set parameters for the ODE solver 
CALL SSET(MXPARM,0.0,PARAM, 1) 
PARAM(1) = HINIT 

PARAM(4) = REAL(MXSTEP) 

PARAM(12) = REAL(METHOD) 


INTIME = MAXTIM / NTIME 
cone eee Initialize queues and Set initial time 


DO 10 1=1,QUANT 
Y(I) = YO(I) 

CONTINUE 

DO 15 I=1,QUANT*QUANT 
Y(QUANT+I) = 0.0 

CONTINUE 

IDO=1 

i INTEGRATE SYSTEM OF ODE 

DO 200 ISTEP =NINT(X0) + INTIME, MAXTIM, INTIME 
XEND = REAL(ISTEP) 


fl 


111 


150 


CALL IVPAG(IDO,QUANT + QUANT*QUANT,FCN,FCNJ,AMAT, 
X,XEND,TOL,PARAM,.Y) 
DO 100 1=1,QUANT 
IF(Y(I).LT.0.0) Y(I) = 0.0 
OK(I) = REAL(K(D) - Y(I)*A 
CONTINUE 


staeiteanas ices FIND COMPONENT WITH LEAST AVAILABILITY 


IMIN = ISMIN(QUANT,OK,1) 


sosceanenetes tae FIND PRESENT AVAILABILITY 


AVNOW = REAL(AC) 
IF(OK(IMIN) .LT. AVNOW) AVNOW = OK(IMIN) 
AVERAG(ISTEP/INTIME) = AVNOW 


een cisenesen COMPUTE STANDARD DEVIATION 


STDDEV(ISTEP /INTIME) = SORT(ABS(A* Y(IPOS(IMIN,IMIN)))) 
WRITE(11,111)ISTEP,(K(I)-Y(1)*A,I = 1,QUANT) 
WRITE(11,111)ISTEP,(SORT(ABS(A*Y(IPOS(LI)))),I = 1,QUANT) 
FORMAT(1X,I3,10(1X,F5.1)) 
DO 150 1=1,QUANT 
DMEAN(ISTEP/INTIME,]) = K(1) - Y(1)*A 
DSTD(ISTEP/INTIME,I) = SQRT(ABS(A* Y(IPOS(LI)))) 
CONTINUE 


200 CONTINUE 


5 


ilcccunptieteeenese RELEASE WORKSPACE 
IDO=3 
CALL IVPAG(IDO,QUANT+ QUANT*QUANT,FCN,FCNJ,AMAT, 
& 
RETURN 
END 


X,XEND,TOL,PARAM,.Y) 


KEKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK 


SUBROUTINE FCN(N,X,Y,YPRIME) 


KKKKKKKKKKKKKEKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK KKK KKK KK 


= 


INTEGER N 
REAL X,Y(N),YPRIME(N) 


INTEGER QUANT 
PARAMETER (QUANT = 10 ) 


Siena MODEL VARIABLES 

INTEGER ONREP,DOWN(QUANT),K(QUANT),NUMMIN,AC 
INTEGER P,AV,TOTQU,A VMIN 

REAL LB(QUANT),NU(QUANT),W(QUANT) 
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COMMON/MODEL1/ONREP,DOWN,K,NUMMIN, AC 
& P,AV,TOTQU,AVMIN 
COMMON/MODEL2/LB,NU,W 


eee DIFFUSION VARIABLES 
REAL MU(QUANT),ALPH(QUANT), Y0(QUANT),ACNORM,X0 
COMMON/APPROX/MU,ALPH, Y0,ACNORM,X0 


Rent ee POLICY VARIABLE 
INTEGER POLTYP 
COMMON/POLICY/POLTYP 


eee... LOCAL VARIABLES 

REAL PJ(QUANT),PSUM,A VNOW,Q(QUANT) 

REAL B(QUANT),H1(QUANT,QUANT),H(QUANT,QUANT) 
REAL SUMHS(QUANT),SUMPMU 

INTEGER IMIN,IPOS,1J,L 


er errs Determine availability (AVNOW) 
IMIN = -1 
AVNOW = ACNORM 
DO 10 I[=1,QUANT 
IF ((ALPH(I)-Y(1)) .LT. AVNOW) THEN 
AVNOW = ALPH(I)-Y(1) 
IMIN = I 
END IF 
CONTINUE 
PSUM=0.0 
DO 25 1=1,QUANT 
eo Ae FCFS CASE 
IF(POLTYP .EQ. 2) THEN 
PJ(I) = Y(I)/MU(1) 
eset otues LAIN CASE 
ELSE IF(POLTYP .EQ. 4) THEN 
PJ(I) = (W()/MU())) / (ALPH(I)-Y(1))**P 
ELSE 


WRITE(*,*)” ERROR IN POLICY TYPE’ 
STOP 

END IF 

PSUM = PSUM + PJ(I) 


73 


25 CONTINUE 
ee COMPUTE ADAPTIVE PRIORITY (Q(i)) 
DO 27 I=1,QUANT 
IF(PSUM .NE. 0.0) THEN 
Q(I) = PJ(I) / PSUM 


ELSE 
Q(1) = 1.0 / REAL(QUANT) 
END IF 
27 CONTINUE 
= _gudeniemnenmcerasenss Compute scaled means 


DO 30 J=1,QUANT 
YPRIME(J)=LB(J) * AVNOW - (MU(J) * Q(J)) 
30 CONTINUE 
eae Compute scaled variance-covariance 
SUMPMU = 0.0 
DO 50 I=1,QUANT 
SUMPMU = SUMPMU + (Q(I)/MU(I)) 
50 CONTINUE 
DO 200 I1=1,QUANT 
BC) = LB(I)*AVNOW + MU(I)*Q(I)*(1.0+2.0* Q(1) 
& *(MU(I)*SUMPMU-1.0)) 
IF(POLTYP .EQ. 2) THEN 
Bee har Eee: FCFS 
IF(Y(1) .NE. 0.0) THEN 
_ NG = (MUG/Y@) * (QM"(.0-M) 


it = (-1.0) * (Q(1)*(1.0-Q(1))) 


END 
ELSE 
oo LAIN CASE 
H1(LI) = (-MU(I)*REAL(P)/(ALPH(1-Y(J))) 
& *(Q(1)*(1.0-Q(D)) 
END IF 


DO 100 J=1,QUANT 
IF(J .NE. I) THEN 
IF(POLTYP .EQ. 2) THEN 
geese FCFS CASE 
IF(Y(I) .NE. 0.0) THEN 
H1(IJ) = (MU(1D)/Y(I))* Q(7)*Q(V) 
ELSE 
H1(LJ) = Q(1) * QQ) 
END IF 
ELSE 
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a ee Ce ee LST AV 
H1(1,J) = (MU(I)*REAL(P)/(ALPH(I)-Y(1))) 
* Q(1) * QV) 
END IF 
END IF 
100 CONTINUE 
200 CONTINUE 
DO 400 I=1,QUANT 2 
SUMHS() = 0.0 
DO 300 J=1,QUANT 
IF(IMIN .GT. 0) .AND. (J .EQ. IMIN)) THEN 
Hd,J) = H1(1,J) - LBC) 
ELSE 
Hd,J) = H1(1,J) 
END IF 
SUMHS(I) = SUMHS() + Hd,J)* YUIPOS(LJ)) 
300 CONTINUE 
400 CONTINUE 
DO 600 I=1,QUANT 
YPRIME(IPOS(,I)) = B(I) + 2.0 * SUMHS(I) 
DO 500 J=1,QUANT 
IF(J .NE. I) THEN 
YPRIME(IPOS(I,J)) = 0.0 
DO 450 L=1,QUANT 
YPRIME(IPOS(I,J)) = YPRIME(IPOS(I,J)) + 
& H(1,L)* YUPOS(J,L)) + H(J,L)* YAPOS(I,L)) 
450 CONTINUE 
END IF 
500 CONTINUE 
600 CONTINUE 
DO 700 I=1,QUANT 
IFCY(I) .LT. 0.0) Y(I) = 0.0 
700 CONTINUE 
RETURN 
END 
INTEGER FUNCTION IPOS(I,J) 
* This function maps matrix-like indices to the vector locations in Y 
INTEGER I,J,QUANT 
PARAMETER (QUANT = 10 ) 
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IPOS = QUANT + I + (J-1) * QUANT 
RETURN 
END 
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APPENDIX D. HEAVY TRAFFIC CONDITION 


In this appendix we show that if the heavy traffic conditions in (1.2) do not hold, 
the diffusion approximation model is profoundly degraded. Consider the sample 


system described in Table 1D. Applying (1.2) to the data yields 


a 
50x) = =»(/8Se< 1 
i=1 VY; 


Table 1D INPUT DATA 
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The heavy traffic condition is obviously violated and we want to verify how this 


situation will influence the accuracy of the diffusion model. Figure 1D displays the 
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expected value of the combat availability for t between 0 and 100. We can see that 


the analytical results show a wide variability when compared to the simulation. 


LAIN Policy 


= 
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Figure 1D. LAIN Policy (Mean) 
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Figure 2D shows that the performance with respect to the calculation of 


variances is extremely poor. 


LAIN Policy 





Figure 2D. LAIN Policy (Standard Deviation) 
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